SIMU=0;CD_REAL_0=numeric()
for(w in 1:10000)
{
y1=PROFILES_G1_0[,w];y2=PROFILES_G2_0[,w]
DIF=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
dif=(x-1)*stap-(y-1)*stap
DIF[x,y]=dif
}}
DIF=sign(DIF)
DEN=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
den=y1[x]*y2[y]*stap^2
DEN[x,y]=den
}}
CD_real=sum(DIF*DEN)/sum(DEN)
CD_REAL_0=c(CD_REAL_0,CD_real)
SIMU=SIMU+1;if(SIMU%%1000==0)print(SIMU)
}
write.table(CD_REAL_0,'CD_REAL_0.txt',row.names=FALSE,col.names=FALSE)

SIMU=0;CD_REAL_1=numeric()
for(w in 1:10000)
{
y1=PROFILES_G1_1[,w];y2=PROFILES_G2_1[,w]
DIF=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
dif=(x-1)*stap-(y-1)*stap
DIF[x,y]=dif
}}
DIF=sign(DIF)
DEN=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
den=y1[x]*y2[y]*stap^2
DEN[x,y]=den
}}
CD_real=sum(DIF*DEN)/sum(DEN)
CD_REAL_1=c(CD_REAL_1,CD_real)
SIMU=SIMU+1;if(SIMU%%1000==0)print(SIMU)
}
write.table(CD_REAL_1,'CD_REAL_1.txt',row.names=FALSE,col.names=FALSE)

SIMU=0;CD_REAL_2=numeric()
for(w in 1:10000)
{
y1=PROFILES_G1_2[,w];y2=PROFILES_G2_2[,w]
DIF=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
dif=(x-1)*stap-(y-1)*stap
DIF[x,y]=dif
}}
DIF=sign(DIF)
DEN=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
den=y1[x]*y2[y]*stap^2
DEN[x,y]=den
}}
CD_real=sum(DIF*DEN)/sum(DEN)
CD_REAL_2=c(CD_REAL_2,CD_real)
SIMU=SIMU+1;if(SIMU%%1000==0)print(SIMU)
}
write.table(CD_REAL_2,'CD_REAL_2.txt',row.names=FALSE,col.names=FALSE)

SIMU=0;CD_REAL_3=numeric()
for(w in 1:10000)
{
y1=PROFILES_G1_3[,w];y2=PROFILES_G2_3[,w]
DIF=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
dif=(x-1)*stap-(y-1)*stap
DIF[x,y]=dif
}}
DIF=sign(DIF)
DEN=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
den=y1[x]*y2[y]*stap^2
DEN[x,y]=den
}}
CD_real=sum(DIF*DEN)/sum(DEN)
CD_REAL_3=c(CD_REAL_3,CD_real)
SIMU=SIMU+1;if(SIMU%%1000==0)print(SIMU)
}
write.table(CD_REAL_3,'CD_REAL_3.txt',row.names=FALSE,col.names=FALSE)

SIMU=0;CD_REAL_4=numeric()
for(w in 1:10000)
{
y1=PROFILES_G1_4[,w];y2=PROFILES_G2_4[,w]
DIF=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
dif=(x-1)*stap-(y-1)*stap
DIF[x,y]=dif
}}
DIF=sign(DIF)
DEN=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
den=y1[x]*y2[y]*stap^2
DEN[x,y]=den
}}
CD_real=sum(DIF*DEN)/sum(DEN)
CD_REAL_4=c(CD_REAL_4,CD_real)
SIMU=SIMU+1;if(SIMU%%1000==0)print(SIMU)
}
write.table(CD_REAL_4,'CD_REAL_4.txt',row.names=FALSE,col.names=FALSE)

SIMU=0;CD_REAL_5=numeric()
for(w in 1:10000)
{
y1=PROFILES_G1_5[,w];y2=PROFILES_G2_5[,w]
DIF=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
dif=(x-1)*stap-(y-1)*stap
DIF[x,y]=dif
}}
DIF=sign(DIF)
DEN=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
den=y1[x]*y2[y]*stap^2
DEN[x,y]=den
}}
CD_real=sum(DIF*DEN)/sum(DEN)
CD_REAL_5=c(CD_REAL_5,CD_real)
SIMU=SIMU+1;if(SIMU%%1000==0)print(SIMU)
}
write.table(CD_REAL_5,'CD_REAL_5.txt',row.names=FALSE,col.names=FALSE)

SIMU=0;CD_REAL_6=numeric()
for(w in 1:10000)
{
y1=PROFILES_G1_6[,w];y2=PROFILES_G2_6[,w]
DIF=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
dif=(x-1)*stap-(y-1)*stap
DIF[x,y]=dif
}}
DIF=sign(DIF)
DEN=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
den=y1[x]*y2[y]*stap^2
DEN[x,y]=den
}}
CD_real=sum(DIF*DEN)/sum(DEN)
CD_REAL_6=c(CD_REAL_6,CD_real)
SIMU=SIMU+1;if(SIMU%%1000==0)print(SIMU)
}
write.table(CD_REAL_6,'CD_REAL_6.txt',row.names=FALSE,col.names=FALSE)

SIMU=0;CD_REAL_7=numeric()
for(w in 1:10000)
{
y1=PROFILES_G1_7[,w];y2=PROFILES_G2_7[,w]
DIF=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
dif=(x-1)*stap-(y-1)*stap
DIF[x,y]=dif
}}
DIF=sign(DIF)
DEN=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
den=y1[x]*y2[y]*stap^2
DEN[x,y]=den
}}
CD_real=sum(DIF*DEN)/sum(DEN)
CD_REAL_7=c(CD_REAL_7,CD_real)
SIMU=SIMU+1;if(SIMU%%1000==0)print(SIMU)
}
write.table(CD_REAL_7,'CD_REAL_7.txt',row.names=FALSE,col.names=FALSE)

SIMU=0;CD_REAL_8=numeric()
for(w in 1:10000)
{
y1=PROFILES_G1_8[,w];y2=PROFILES_G2_8[,w]
DIF=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
dif=(x-1)*stap-(y-1)*stap
DIF[x,y]=dif
}}
DIF=sign(DIF)
DEN=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
den=y1[x]*y2[y]*stap^2
DEN[x,y]=den
}}
CD_real=sum(DIF*DEN)/sum(DEN)
CD_REAL_8=c(CD_REAL_8,CD_real)
SIMU=SIMU+1;if(SIMU%%1000==0)print(SIMU)
}
write.table(CD_REAL_8,'CD_REAL_8.txt',row.names=FALSE,col.names=FALSE)

SIMU=0;CD_REAL_9=numeric()
for(w in 1:10000)
{
y1=PROFILES_G1_9[,w];y2=PROFILES_G2_9[,w]
DIF=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
dif=(x-1)*stap-(y-1)*stap
DIF[x,y]=dif
}}
DIF=sign(DIF)
DEN=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
den=y1[x]*y2[y]*stap^2
DEN[x,y]=den
}}
CD_real=sum(DIF*DEN)/sum(DEN)
CD_REAL_9=c(CD_REAL_9,CD_real)
SIMU=SIMU+1;if(SIMU%%1000==0)print(SIMU)
}
write.table(CD_REAL_9,'CD_REAL_9.txt',row.names=FALSE,col.names=FALSE)

SIMU=0;CD_REAL_10=numeric()
for(w in 1:10000)
{
y1=PROFILES_G1_10[,w];y2=PROFILES_G2_10[,w]
DIF=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
dif=(x-1)*stap-(y-1)*stap
DIF[x,y]=dif
}}
DIF=sign(DIF)
DEN=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
den=y1[x]*y2[y]*stap^2
DEN[x,y]=den
}}
CD_real=sum(DIF*DEN)/sum(DEN)
CD_REAL_10=c(CD_REAL_10,CD_real)
SIMU=SIMU+1;if(SIMU%%1000==0)print(SIMU)
}
write.table(CD_REAL_10,'CD_REAL_10.txt',row.names=FALSE,col.names=FALSE)

SIMU=0;CD_REAL_15=numeric()
for(w in 1:10000)
{
y1=PROFILES_G1_15[,w];y2=PROFILES_G2_15[,w]
DIF=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
dif=(x-1)*stap-(y-1)*stap
DIF[x,y]=dif
}}
DIF=sign(DIF)
DEN=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
den=y1[x]*y2[y]*stap^2
DEN[x,y]=den
}}
CD_real=sum(DIF*DEN)/sum(DEN)
CD_REAL_15=c(CD_REAL_15,CD_real)
SIMU=SIMU+1;if(SIMU%%1000==0)print(SIMU)
}
write.table(CD_REAL_15,'CD_REAL_15.txt',row.names=FALSE,col.names=FALSE)

SIMU=0;CD_REAL_20=numeric()
for(w in 1:10000)
{
y1=PROFILES_G1_20[,w];y2=PROFILES_G2_20[,w]
DIF=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
dif=(x-1)*stap-(y-1)*stap
DIF[x,y]=dif
}}
DIF=sign(DIF)
DEN=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
den=y1[x]*y2[y]*stap^2
DEN[x,y]=den
}}
CD_real=sum(DIF*DEN)/sum(DEN)
CD_REAL_20=c(CD_REAL_20,CD_real)
SIMU=SIMU+1;if(SIMU%%1000==0)print(SIMU)
}
write.table(CD_REAL_20,'CD_REAL_20.txt',row.names=FALSE,col.names=FALSE)

SIMU=0;CD_REAL_25=numeric()
for(w in 1:10000)
{
y1=PROFILES_G1_25[,w];y2=PROFILES_G2_25[,w]
DIF=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
dif=(x-1)*stap-(y-1)*stap
DIF[x,y]=dif
}}
DIF=sign(DIF)
DEN=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
den=y1[x]*y2[y]*stap^2
DEN[x,y]=den
}}
CD_real=sum(DIF*DEN)/sum(DEN)
CD_REAL_25=c(CD_REAL_25,CD_real)
SIMU=SIMU+1;if(SIMU%%1000==0)print(SIMU)
}
write.table(CD_REAL_25,'CD_REAL_25.txt',row.names=FALSE,col.names=FALSE)

SIMU=0;CD_REAL_30=numeric()
for(w in 1:10000)
{
y1=PROFILES_G1_30[,w];y2=PROFILES_G2_30[,w]
DIF=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
dif=(x-1)*stap-(y-1)*stap
DIF[x,y]=dif
}}
DIF=sign(DIF)
DEN=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
den=y1[x]*y2[y]*stap^2
DEN[x,y]=den
}}
CD_real=sum(DIF*DEN)/sum(DEN)
CD_REAL_30=c(CD_REAL_30,CD_real)
SIMU=SIMU+1;if(SIMU%%1000==0)print(SIMU)
}
write.table(CD_REAL_30,'CD_REAL_30.txt',row.names=FALSE,col.names=FALSE)

SIMU=0;CD_REAL_35=numeric()
for(w in 1:10000)
{
y1=PROFILES_G1_35[,w];y2=PROFILES_G2_35[,w]
DIF=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
dif=(x-1)*stap-(y-1)*stap
DIF[x,y]=dif
}}
DIF=sign(DIF)
DEN=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
den=y1[x]*y2[y]*stap^2
DEN[x,y]=den
}}
CD_real=sum(DIF*DEN)/sum(DEN)
CD_REAL_35=c(CD_REAL_35,CD_real)
SIMU=SIMU+1;if(SIMU%%1000==0)print(SIMU)
}
write.table(CD_REAL_35,'CD_REAL_35.txt',row.names=FALSE,col.names=FALSE)

SIMU=0;CD_REAL_40=numeric()
for(w in 1:10000)
{
y1=PROFILES_G1_40[,w];y2=PROFILES_G2_40[,w]
DIF=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
dif=(x-1)*stap-(y-1)*stap
DIF[x,y]=dif
}}
DIF=sign(DIF)
DEN=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
den=y1[x]*y2[y]*stap^2
DEN[x,y]=den
}}
CD_real=sum(DIF*DEN)/sum(DEN)
CD_REAL_40=c(CD_REAL_40,CD_real)
SIMU=SIMU+1;if(SIMU%%1000==0)print(SIMU)
}
write.table(CD_REAL_40,'CD_REAL_40.txt',row.names=FALSE,col.names=FALSE)

SIMU=0;CD_REAL_45=numeric()
for(w in 1:10000)
{
y1=PROFILES_G1_45[,w];y2=PROFILES_G2_45[,w]
DIF=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
dif=(x-1)*stap-(y-1)*stap
DIF[x,y]=dif
}}
DIF=sign(DIF)
DEN=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
den=y1[x]*y2[y]*stap^2
DEN[x,y]=den
}}
CD_real=sum(DIF*DEN)/sum(DEN)
CD_REAL_45=c(CD_REAL_45,CD_real)
SIMU=SIMU+1;if(SIMU%%1000==0)print(SIMU)
}
write.table(CD_REAL_45,'CD_REAL_45.txt',row.names=FALSE,col.names=FALSE)

SIMU=0;CD_REAL_50=numeric()
for(w in 1:10000)
{
y1=PROFILES_G1_50[,w];y2=PROFILES_G2_50[,w]
DIF=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
dif=(x-1)*stap-(y-1)*stap
DIF[x,y]=dif
}}
DIF=sign(DIF)
DEN=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
den=y1[x]*y2[y]*stap^2
DEN[x,y]=den
}}
CD_real=sum(DIF*DEN)/sum(DEN)
CD_REAL_50=c(CD_REAL_50,CD_real)
SIMU=SIMU+1;if(SIMU%%1000==0)print(SIMU)
}
write.table(CD_REAL_50,'CD_REAL_50.txt',row.names=FALSE,col.names=FALSE)

SIMU=0;CD_REAL_55=numeric()
for(w in 1:10000)
{
y1=PROFILES_G1_55[,w];y2=PROFILES_G2_55[,w]
DIF=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
dif=(x-1)*stap-(y-1)*stap
DIF[x,y]=dif
}}
DIF=sign(DIF)
DEN=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
den=y1[x]*y2[y]*stap^2
DEN[x,y]=den
}}
CD_real=sum(DIF*DEN)/sum(DEN)
CD_REAL_55=c(CD_REAL_55,CD_real)
SIMU=SIMU+1;if(SIMU%%1000==0)print(SIMU)
}
write.table(CD_REAL_55,'CD_REAL_55.txt',row.names=FALSE,col.names=FALSE)

SIMU=0;CD_REAL_60=numeric()
for(w in 1:10000)
{
y1=PROFILES_G1_60[,w];y2=PROFILES_G2_60[,w]
DIF=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
dif=(x-1)*stap-(y-1)*stap
DIF[x,y]=dif
}}
DIF=sign(DIF)
DEN=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
den=y1[x]*y2[y]*stap^2
DEN[x,y]=den
}}
CD_real=sum(DIF*DEN)/sum(DEN)
CD_REAL_60=c(CD_REAL_60,CD_real)
SIMU=SIMU+1;if(SIMU%%1000==0)print(SIMU)
}
write.table(CD_REAL_60,'CD_REAL_60.txt',row.names=FALSE,col.names=FALSE)

SIMU=0;CD_REAL_65=numeric()
for(w in 1:10000)
{
y1=PROFILES_G1_65[,w];y2=PROFILES_G2_65[,w]
DIF=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
dif=(x-1)*stap-(y-1)*stap
DIF[x,y]=dif
}}
DIF=sign(DIF)
DEN=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
den=y1[x]*y2[y]*stap^2
DEN[x,y]=den
}}
CD_real=sum(DIF*DEN)/sum(DEN)
CD_REAL_65=c(CD_REAL_65,CD_real)
SIMU=SIMU+1;if(SIMU%%1000==0)print(SIMU)
}
write.table(CD_REAL_65,'CD_REAL_65.txt',row.names=FALSE,col.names=FALSE)

SIMU=0;CD_REAL_70=numeric()
for(w in 1:10000)
{
y1=PROFILES_G1_70[,w];y2=PROFILES_G2_70[,w]
DIF=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
dif=(x-1)*stap-(y-1)*stap
DIF[x,y]=dif
}}
DIF=sign(DIF)
DEN=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
den=y1[x]*y2[y]*stap^2
DEN[x,y]=den
}}
CD_real=sum(DIF*DEN)/sum(DEN)
CD_REAL_70=c(CD_REAL_70,CD_real)
SIMU=SIMU+1;if(SIMU%%1000==0)print(SIMU)
}
write.table(CD_REAL_70,'CD_REAL_70.txt',row.names=FALSE,col.names=FALSE)

SIMU=0;CD_REAL_75=numeric()
for(w in 1:10000)
{
y1=PROFILES_G1_75[,w];y2=PROFILES_G2_75[,w]
DIF=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
dif=(x-1)*stap-(y-1)*stap
DIF[x,y]=dif
}}
DIF=sign(DIF)
DEN=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
den=y1[x]*y2[y]*stap^2
DEN[x,y]=den
}}
CD_real=sum(DIF*DEN)/sum(DEN)
CD_REAL_75=c(CD_REAL_75,CD_real)
SIMU=SIMU+1;if(SIMU%%1000==0)print(SIMU)
}
write.table(CD_REAL_75,'CD_REAL_75.txt',row.names=FALSE,col.names=FALSE)

SIMU=0;CD_REAL_80=numeric()
for(w in 1:10000)
{
y1=PROFILES_G1_80[,w];y2=PROFILES_G2_80[,w]
DIF=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
dif=(x-1)*stap-(y-1)*stap
DIF[x,y]=dif
}}
DIF=sign(DIF)
DEN=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
den=y1[x]*y2[y]*stap^2
DEN[x,y]=den
}}
CD_real=sum(DIF*DEN)/sum(DEN)
CD_REAL_80=c(CD_REAL_80,CD_real)
SIMU=SIMU+1;if(SIMU%%1000==0)print(SIMU)
}
write.table(CD_REAL_80,'CD_REAL_80.txt',row.names=FALSE,col.names=FALSE)

SIMU=0;CD_REAL_85=numeric()
for(w in 1:10000)
{
y1=PROFILES_G1_85[,w];y2=PROFILES_G2_85[,w]
DIF=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
dif=(x-1)*stap-(y-1)*stap
DIF[x,y]=dif
}}
DIF=sign(DIF)
DEN=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
den=y1[x]*y2[y]*stap^2
DEN[x,y]=den
}}
CD_real=sum(DIF*DEN)/sum(DEN)
CD_REAL_85=c(CD_REAL_85,CD_real)
SIMU=SIMU+1;if(SIMU%%1000==0)print(SIMU)
}
write.table(CD_REAL_85,'CD_REAL_85.txt',row.names=FALSE,col.names=FALSE)

SIMU=0;CD_REAL_90=numeric()
for(w in 1:10000)
{
y1=PROFILES_G1_90[,w];y2=PROFILES_G2_90[,w]
DIF=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
dif=(x-1)*stap-(y-1)*stap
DIF[x,y]=dif
}}
DIF=sign(DIF)
DEN=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
den=y1[x]*y2[y]*stap^2
DEN[x,y]=den
}}
CD_real=sum(DIF*DEN)/sum(DEN)
CD_REAL_90=c(CD_REAL_90,CD_real)
SIMU=SIMU+1;if(SIMU%%1000==0)print(SIMU)
}
write.table(CD_REAL_90,'CD_REAL_90.txt',row.names=FALSE,col.names=FALSE)

SIMU=0;CD_REAL_95=numeric()
for(w in 1:10000)
{
y1=PROFILES_G1_95[,w];y2=PROFILES_G2_95[,w]
DIF=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
dif=(x-1)*stap-(y-1)*stap
DIF[x,y]=dif
}}
DIF=sign(DIF)
DEN=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
den=y1[x]*y2[y]*stap^2
DEN[x,y]=den
}}
CD_real=sum(DIF*DEN)/sum(DEN)
CD_REAL_95=c(CD_REAL_95,CD_real)
SIMU=SIMU+1;if(SIMU%%1000==0)print(SIMU)
}
write.table(CD_REAL_95,'CD_REAL_95.txt',row.names=FALSE,col.names=FALSE)

SIMU=0;CD_REAL_100=numeric()
for(w in 1:10000)
{
y1=PROFILES_G1_100[,w];y2=PROFILES_G2_100[,w]
DIF=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
dif=(x-1)*stap-(y-1)*stap
DIF[x,y]=dif
}}
DIF=sign(DIF)
DEN=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
den=y1[x]*y2[y]*stap^2
DEN[x,y]=den
}}
CD_real=sum(DIF*DEN)/sum(DEN)
CD_REAL_100=c(CD_REAL_100,CD_real)
SIMU=SIMU+1;if(SIMU%%1000==0)print(SIMU)
}
write.table(CD_REAL_100,'CD_REAL_100.txt',row.names=FALSE,col.names=FALSE)

SIMU=0;CD_REAL_110=numeric()
for(w in 1:10000)
{
y1=PROFILES_G1_110[,w];y2=PROFILES_G2_110[,w]
DIF=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
dif=(x-1)*stap-(y-1)*stap
DIF[x,y]=dif
}}
DIF=sign(DIF)
DEN=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
den=y1[x]*y2[y]*stap^2
DEN[x,y]=den
}}
CD_real=sum(DIF*DEN)/sum(DEN)
CD_REAL_110=c(CD_REAL_110,CD_real)
SIMU=SIMU+1;if(SIMU%%1000==0)print(SIMU)
}
write.table(CD_REAL_110,'CD_REAL_110.txt',row.names=FALSE,col.names=FALSE)

SIMU=0;CD_REAL_120=numeric()
for(w in 1:10000)
{
y1=PROFILES_G1_120[,w];y2=PROFILES_G2_120[,w]
DIF=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
dif=(x-1)*stap-(y-1)*stap
DIF[x,y]=dif
}}
DIF=sign(DIF)
DEN=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
den=y1[x]*y2[y]*stap^2
DEN[x,y]=den
}}
CD_real=sum(DIF*DEN)/sum(DEN)
CD_REAL_120=c(CD_REAL_120,CD_real)
SIMU=SIMU+1;if(SIMU%%1000==0)print(SIMU)
}
write.table(CD_REAL_120,'CD_REAL_120.txt',row.names=FALSE,col.names=FALSE)

SIMU=0;CD_REAL_130=numeric()
for(w in 1:10000)
{
y1=PROFILES_G1_130[,w];y2=PROFILES_G2_130[,w]
DIF=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
dif=(x-1)*stap-(y-1)*stap
DIF[x,y]=dif
}}
DIF=sign(DIF)
DEN=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
den=y1[x]*y2[y]*stap^2
DEN[x,y]=den
}}
CD_real=sum(DIF*DEN)/sum(DEN)
CD_REAL_130=c(CD_REAL_130,CD_real)
SIMU=SIMU+1;if(SIMU%%1000==0)print(SIMU)
}
write.table(CD_REAL_130,'CD_REAL_130.txt',row.names=FALSE,col.names=FALSE)

SIMU=0;CD_REAL_140=numeric()
for(w in 1:10000)
{
y1=PROFILES_G1_140[,w];y2=PROFILES_G2_140[,w]
DIF=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
dif=(x-1)*stap-(y-1)*stap
DIF[x,y]=dif
}}
DIF=sign(DIF)
DEN=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
den=y1[x]*y2[y]*stap^2
DEN[x,y]=den
}}
CD_real=sum(DIF*DEN)/sum(DEN)
CD_REAL_140=c(CD_REAL_140,CD_real)
SIMU=SIMU+1;if(SIMU%%1000==0)print(SIMU)
}
write.table(CD_REAL_140,'CD_REAL_140.txt',row.names=FALSE,col.names=FALSE)

SIMU=0;CD_REAL_150=numeric()
for(w in 1:10000)
{
y1=PROFILES_G1_150[,w];y2=PROFILES_G2_150[,w]
DIF=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
dif=(x-1)*stap-(y-1)*stap
DIF[x,y]=dif
}}
DIF=sign(DIF)
DEN=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
den=y1[x]*y2[y]*stap^2
DEN[x,y]=den
}}
CD_real=sum(DIF*DEN)/sum(DEN)
CD_REAL_150=c(CD_REAL_150,CD_real)
SIMU=SIMU+1;if(SIMU%%1000==0)print(SIMU)
}
write.table(CD_REAL_150,'CD_REAL_150.txt',row.names=FALSE,col.names=FALSE)

SIMU=0;CD_REAL_160=numeric()
for(w in 1:10000)
{
y1=PROFILES_G1_160[,w];y2=PROFILES_G2_160[,w]
DIF=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
dif=(x-1)*stap-(y-1)*stap
DIF[x,y]=dif
}}
DIF=sign(DIF)
DEN=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
den=y1[x]*y2[y]*stap^2
DEN[x,y]=den
}}
CD_real=sum(DIF*DEN)/sum(DEN)
CD_REAL_160=c(CD_REAL_160,CD_real)
SIMU=SIMU+1;if(SIMU%%1000==0)print(SIMU)
}
write.table(CD_REAL_160,'CD_REAL_160.txt',row.names=FALSE,col.names=FALSE)

SIMU=0;CD_REAL_170=numeric()
for(w in 1:10000)
{
y1=PROFILES_G1_170[,w];y2=PROFILES_G2_170[,w]
DIF=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
dif=(x-1)*stap-(y-1)*stap
DIF[x,y]=dif
}}
DIF=sign(DIF)
DEN=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
den=y1[x]*y2[y]*stap^2
DEN[x,y]=den
}}
CD_real=sum(DIF*DEN)/sum(DEN)
CD_REAL_170=c(CD_REAL_170,CD_real)
SIMU=SIMU+1;if(SIMU%%1000==0)print(SIMU)
}
write.table(CD_REAL_170,'CD_REAL_170.txt',row.names=FALSE,col.names=FALSE)

SIMU=0;CD_REAL_180=numeric()
for(w in 1:10000)
{
y1=PROFILES_G1_180[,w];y2=PROFILES_G2_180[,w]
DIF=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
dif=(x-1)*stap-(y-1)*stap
DIF[x,y]=dif
}}
DIF=sign(DIF)
DEN=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
den=y1[x]*y2[y]*stap^2
DEN[x,y]=den
}}
CD_real=sum(DIF*DEN)/sum(DEN)
CD_REAL_180=c(CD_REAL_180,CD_real)
SIMU=SIMU+1;if(SIMU%%1000==0)print(SIMU)
}
write.table(CD_REAL_180,'CD_REAL_180.txt',row.names=FALSE,col.names=FALSE)

SIMU=0;CD_REAL_190=numeric()
for(w in 1:10000)
{
y1=PROFILES_G1_190[,w];y2=PROFILES_G2_190[,w]
DIF=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
dif=(x-1)*stap-(y-1)*stap
DIF[x,y]=dif
}}
DIF=sign(DIF)
DEN=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
den=y1[x]*y2[y]*stap^2
DEN[x,y]=den
}}
CD_real=sum(DIF*DEN)/sum(DEN)
CD_REAL_190=c(CD_REAL_190,CD_real)
SIMU=SIMU+1;if(SIMU%%1000==0)print(SIMU)
}
write.table(CD_REAL_190,'CD_REAL_190.txt',row.names=FALSE,col.names=FALSE)

SIMU=0;CD_REAL_200=numeric()
for(w in 1:10000)
{
y1=PROFILES_G1_200[,w];y2=PROFILES_G2_200[,w]
DIF=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
dif=(x-1)*stap-(y-1)*stap
DIF[x,y]=dif
}}
DIF=sign(DIF)
DEN=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
den=y1[x]*y2[y]*stap^2
DEN[x,y]=den
}}
CD_real=sum(DIF*DEN)/sum(DEN)
CD_REAL_200=c(CD_REAL_200,CD_real)
SIMU=SIMU+1;if(SIMU%%1000==0)print(SIMU)
}
write.table(CD_REAL_200,'CD_REAL_200.txt',row.names=FALSE,col.names=FALSE)

SIMU=0;CD_REAL_210=numeric()
for(w in 1:10000)
{
y1=PROFILES_G1_210[,w];y2=PROFILES_G2_210[,w]
DIF=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
dif=(x-1)*stap-(y-1)*stap
DIF[x,y]=dif
}}
DIF=sign(DIF)
DEN=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
den=y1[x]*y2[y]*stap^2
DEN[x,y]=den
}}
CD_real=sum(DIF*DEN)/sum(DEN)
CD_REAL_210=c(CD_REAL_210,CD_real)
SIMU=SIMU+1;if(SIMU%%1000==0)print(SIMU)
}
write.table(CD_REAL_210,'CD_REAL_210.txt',row.names=FALSE,col.names=FALSE)

SIMU=0;CD_REAL_220=numeric()
for(w in 1:10000)
{
y1=PROFILES_G1_220[,w];y2=PROFILES_G2_220[,w]
DIF=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
dif=(x-1)*stap-(y-1)*stap
DIF[x,y]=dif
}}
DIF=sign(DIF)
DEN=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
den=y1[x]*y2[y]*stap^2
DEN[x,y]=den
}}
CD_real=sum(DIF*DEN)/sum(DEN)
CD_REAL_220=c(CD_REAL_220,CD_real)
SIMU=SIMU+1;if(SIMU%%1000==0)print(SIMU)
}
write.table(CD_REAL_220,'CD_REAL_220.txt',row.names=FALSE,col.names=FALSE)

SIMU=0;CD_REAL_230=numeric()
for(w in 1:10000)
{
y1=PROFILES_G1_230[,w];y2=PROFILES_G2_230[,w]
DIF=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
dif=(x-1)*stap-(y-1)*stap
DIF[x,y]=dif
}}
DIF=sign(DIF)
DEN=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
den=y1[x]*y2[y]*stap^2
DEN[x,y]=den
}}
CD_real=sum(DIF*DEN)/sum(DEN)
CD_REAL_230=c(CD_REAL_230,CD_real)
SIMU=SIMU+1;if(SIMU%%1000==0)print(SIMU)
}
write.table(CD_REAL_230,'CD_REAL_230.txt',row.names=FALSE,col.names=FALSE)

SIMU=0;CD_REAL_240=numeric()
for(w in 1:10000)
{
y1=PROFILES_G1_240[,w];y2=PROFILES_G2_240[,w]
DIF=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
dif=(x-1)*stap-(y-1)*stap
DIF[x,y]=dif
}}
DIF=sign(DIF)
DEN=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
den=y1[x]*y2[y]*stap^2
DEN[x,y]=den
}}
CD_real=sum(DIF*DEN)/sum(DEN)
CD_REAL_240=c(CD_REAL_240,CD_real)
SIMU=SIMU+1;if(SIMU%%1000==0)print(SIMU)
}
write.table(CD_REAL_240,'CD_REAL_240.txt',row.names=FALSE,col.names=FALSE)

SIMU=0;CD_REAL_250=numeric()
for(w in 1:10000)
{
y1=PROFILES_G1_250[,w];y2=PROFILES_G2_250[,w]
DIF=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
dif=(x-1)*stap-(y-1)*stap
DIF[x,y]=dif
}}
DIF=sign(DIF)
DEN=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
den=y1[x]*y2[y]*stap^2
DEN[x,y]=den
}}
CD_real=sum(DIF*DEN)/sum(DEN)
CD_REAL_250=c(CD_REAL_250,CD_real)
SIMU=SIMU+1;if(SIMU%%1000==0)print(SIMU)
}
write.table(CD_REAL_250,'CD_REAL_250.txt',row.names=FALSE,col.names=FALSE)

SIMU=0;CD_REAL_260=numeric()
for(w in 1:10000)
{
y1=PROFILES_G1_260[,w];y2=PROFILES_G2_260[,w]
DIF=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
dif=(x-1)*stap-(y-1)*stap
DIF[x,y]=dif
}}
DIF=sign(DIF)
DEN=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
den=y1[x]*y2[y]*stap^2
DEN[x,y]=den
}}
CD_real=sum(DIF*DEN)/sum(DEN)
CD_REAL_260=c(CD_REAL_260,CD_real)
SIMU=SIMU+1;if(SIMU%%1000==0)print(SIMU)
}
write.table(CD_REAL_260,'CD_REAL_260.txt',row.names=FALSE,col.names=FALSE)

SIMU=0;CD_REAL_270=numeric()
for(w in 1:10000)
{
y1=PROFILES_G1_270[,w];y2=PROFILES_G2_270[,w]
DIF=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
dif=(x-1)*stap-(y-1)*stap
DIF[x,y]=dif
}}
DIF=sign(DIF)
DEN=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
den=y1[x]*y2[y]*stap^2
DEN[x,y]=den
}}
CD_real=sum(DIF*DEN)/sum(DEN)
CD_REAL_270=c(CD_REAL_270,CD_real)
SIMU=SIMU+1;if(SIMU%%1000==0)print(SIMU)
}
write.table(CD_REAL_270,'CD_REAL_270.txt',row.names=FALSE,col.names=FALSE)

SIMU=0;CD_REAL_280=numeric()
for(w in 1:10000)
{
y1=PROFILES_G1_280[,w];y2=PROFILES_G2_280[,w]
DIF=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
dif=(x-1)*stap-(y-1)*stap
DIF[x,y]=dif
}}
DIF=sign(DIF)
DEN=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
den=y1[x]*y2[y]*stap^2
DEN[x,y]=den
}}
CD_real=sum(DIF*DEN)/sum(DEN)
CD_REAL_280=c(CD_REAL_280,CD_real)
SIMU=SIMU+1;if(SIMU%%1000==0)print(SIMU)
}
write.table(CD_REAL_280,'CD_REAL_280.txt',row.names=FALSE,col.names=FALSE)

SIMU=0;CD_REAL_290=numeric()
for(w in 1:10000)
{
y1=PROFILES_G1_290[,w];y2=PROFILES_G2_290[,w]
DIF=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
dif=(x-1)*stap-(y-1)*stap
DIF[x,y]=dif
}}
DIF=sign(DIF)
DEN=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
den=y1[x]*y2[y]*stap^2
DEN[x,y]=den
}}
CD_real=sum(DIF*DEN)/sum(DEN)
CD_REAL_290=c(CD_REAL_290,CD_real)
SIMU=SIMU+1;if(SIMU%%1000==0)print(SIMU)
}
write.table(CD_REAL_290,'CD_REAL_290.txt',row.names=FALSE,col.names=FALSE)

SIMU=0;CD_REAL_300=numeric()
for(w in 1:10000)
{
y1=PROFILES_G1_300[,w];y2=PROFILES_G2_300[,w]
DIF=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
dif=(x-1)*stap-(y-1)*stap
DIF[x,y]=dif
}}
DIF=sign(DIF)
DEN=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
den=y1[x]*y2[y]*stap^2
DEN[x,y]=den
}}
CD_real=sum(DIF*DEN)/sum(DEN)
CD_REAL_300=c(CD_REAL_300,CD_real)
SIMU=SIMU+1;if(SIMU%%1000==0)print(SIMU)
}
write.table(CD_REAL_300,'CD_REAL_300.txt',row.names=FALSE,col.names=FALSE)

SIMU=0;CD_REAL_310=numeric()
for(w in 1:10000)
{
y1=PROFILES_G1_310[,w];y2=PROFILES_G2_310[,w]
DIF=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
dif=(x-1)*stap-(y-1)*stap
DIF[x,y]=dif
}}
DIF=sign(DIF)
DEN=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
den=y1[x]*y2[y]*stap^2
DEN[x,y]=den
}}
CD_real=sum(DIF*DEN)/sum(DEN)
CD_REAL_310=c(CD_REAL_310,CD_real)
SIMU=SIMU+1;if(SIMU%%1000==0)print(SIMU)
}
write.table(CD_REAL_310,'CD_REAL_310.txt',row.names=FALSE,col.names=FALSE)

SIMU=0;CD_REAL_320=numeric()
for(w in 1:10000)
{
y1=PROFILES_G1_320[,w];y2=PROFILES_G2_320[,w]
DIF=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
dif=(x-1)*stap-(y-1)*stap
DIF[x,y]=dif
}}
DIF=sign(DIF)
DEN=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
den=y1[x]*y2[y]*stap^2
DEN[x,y]=den
}}
CD_real=sum(DIF*DEN)/sum(DEN)
CD_REAL_320=c(CD_REAL_320,CD_real)
SIMU=SIMU+1;if(SIMU%%1000==0)print(SIMU)
}
write.table(CD_REAL_320,'CD_REAL_320.txt',row.names=FALSE,col.names=FALSE)

SIMU=0;CD_REAL_330=numeric()
for(w in 1:10000)
{
y1=PROFILES_G1_330[,w];y2=PROFILES_G2_330[,w]
DIF=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
dif=(x-1)*stap-(y-1)*stap
DIF[x,y]=dif
}}
DIF=sign(DIF)
DEN=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
den=y1[x]*y2[y]*stap^2
DEN[x,y]=den
}}
CD_real=sum(DIF*DEN)/sum(DEN)
CD_REAL_330=c(CD_REAL_330,CD_real)
SIMU=SIMU+1;if(SIMU%%1000==0)print(SIMU)
}
write.table(CD_REAL_330,'CD_REAL_330.txt',row.names=FALSE,col.names=FALSE)

SIMU=0;CD_REAL_340=numeric()
for(w in 1:10000)
{
y1=PROFILES_G1_340[,w];y2=PROFILES_G2_340[,w]
DIF=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
dif=(x-1)*stap-(y-1)*stap
DIF[x,y]=dif
}}
DIF=sign(DIF)
DEN=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
den=y1[x]*y2[y]*stap^2
DEN[x,y]=den
}}
CD_real=sum(DIF*DEN)/sum(DEN)
CD_REAL_340=c(CD_REAL_340,CD_real)
SIMU=SIMU+1;if(SIMU%%1000==0)print(SIMU)
}
write.table(CD_REAL_340,'CD_REAL_340.txt',row.names=FALSE,col.names=FALSE)

SIMU=0;CD_REAL_350=numeric()
for(w in 1:10000)
{
y1=PROFILES_G1_350[,w];y2=PROFILES_G2_350[,w]
DIF=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
dif=(x-1)*stap-(y-1)*stap
DIF[x,y]=dif
}}
DIF=sign(DIF)
DEN=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
den=y1[x]*y2[y]*stap^2
DEN[x,y]=den
}}
CD_real=sum(DIF*DEN)/sum(DEN)
CD_REAL_350=c(CD_REAL_350,CD_real)
SIMU=SIMU+1;if(SIMU%%1000==0)print(SIMU)
}
write.table(CD_REAL_350,'CD_REAL_350.txt',row.names=FALSE,col.names=FALSE)

SIMU=0;CD_REAL_360=numeric()
for(w in 1:10000)
{
y1=PROFILES_G1_360[,w];y2=PROFILES_G2_360[,w]
DIF=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
dif=(x-1)*stap-(y-1)*stap
DIF[x,y]=dif
}}
DIF=sign(DIF)
DEN=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
den=y1[x]*y2[y]*stap^2
DEN[x,y]=den
}}
CD_real=sum(DIF*DEN)/sum(DEN)
CD_REAL_360=c(CD_REAL_360,CD_real)
SIMU=SIMU+1;if(SIMU%%1000==0)print(SIMU)
}
write.table(CD_REAL_360,'CD_REAL_360.txt',row.names=FALSE,col.names=FALSE)

SIMU=0;CD_REAL_370=numeric()
for(w in 1:10000)
{
y1=PROFILES_G1_370[,w];y2=PROFILES_G2_370[,w]
DIF=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
dif=(x-1)*stap-(y-1)*stap
DIF[x,y]=dif
}}
DIF=sign(DIF)
DEN=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
den=y1[x]*y2[y]*stap^2
DEN[x,y]=den
}}
CD_real=sum(DIF*DEN)/sum(DEN)
CD_REAL_370=c(CD_REAL_370,CD_real)
SIMU=SIMU+1;if(SIMU%%1000==0)print(SIMU)
}
write.table(CD_REAL_370,'CD_REAL_370.txt',row.names=FALSE,col.names=FALSE)

SIMU=0;CD_REAL_380=numeric()
for(w in 1:10000)
{
y1=PROFILES_G1_380[,w];y2=PROFILES_G2_380[,w]
DIF=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
dif=(x-1)*stap-(y-1)*stap
DIF[x,y]=dif
}}
DIF=sign(DIF)
DEN=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
den=y1[x]*y2[y]*stap^2
DEN[x,y]=den
}}
CD_real=sum(DIF*DEN)/sum(DEN)
CD_REAL_380=c(CD_REAL_380,CD_real)
SIMU=SIMU+1;if(SIMU%%1000==0)print(SIMU)
}
write.table(CD_REAL_380,'CD_REAL_380.txt',row.names=FALSE,col.names=FALSE)

SIMU=0;CD_REAL_390=numeric()
for(w in 1:10000)
{
y1=PROFILES_G1_390[,w];y2=PROFILES_G2_390[,w]
DIF=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
dif=(x-1)*stap-(y-1)*stap
DIF[x,y]=dif
}}
DIF=sign(DIF)
DEN=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
den=y1[x]*y2[y]*stap^2
DEN[x,y]=den
}}
CD_real=sum(DIF*DEN)/sum(DEN)
CD_REAL_390=c(CD_REAL_390,CD_real)
SIMU=SIMU+1;if(SIMU%%1000==0)print(SIMU)
}
write.table(CD_REAL_390,'CD_REAL_390.txt',row.names=FALSE,col.names=FALSE)

SIMU=0;CD_REAL_400=numeric()
for(w in 1:10000)
{
y1=PROFILES_G1_400[,w];y2=PROFILES_G2_400[,w]
DIF=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
dif=(x-1)*stap-(y-1)*stap
DIF[x,y]=dif
}}
DIF=sign(DIF)
DEN=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
den=y1[x]*y2[y]*stap^2
DEN[x,y]=den
}}
CD_real=sum(DIF*DEN)/sum(DEN)
CD_REAL_400=c(CD_REAL_400,CD_real)
SIMU=SIMU+1;if(SIMU%%1000==0)print(SIMU)
}
write.table(CD_REAL_400,'CD_REAL_400.txt',row.names=FALSE,col.names=FALSE)

SIMU=0;CD_REAL_410=numeric()
for(w in 1:10000)
{
y1=PROFILES_G1_410[,w];y2=PROFILES_G2_410[,w]
DIF=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
dif=(x-1)*stap-(y-1)*stap
DIF[x,y]=dif
}}
DIF=sign(DIF)
DEN=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
den=y1[x]*y2[y]*stap^2
DEN[x,y]=den
}}
CD_real=sum(DIF*DEN)/sum(DEN)
CD_REAL_410=c(CD_REAL_410,CD_real)
SIMU=SIMU+1;if(SIMU%%1000==0)print(SIMU)
}
write.table(CD_REAL_410,'CD_REAL_410.txt',row.names=FALSE,col.names=FALSE)

SIMU=0;CD_REAL_420=numeric()
for(w in 1:10000)
{
y1=PROFILES_G1_420[,w];y2=PROFILES_G2_420[,w]
DIF=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
dif=(x-1)*stap-(y-1)*stap
DIF[x,y]=dif
}}
DIF=sign(DIF)
DEN=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
den=y1[x]*y2[y]*stap^2
DEN[x,y]=den
}}
CD_real=sum(DIF*DEN)/sum(DEN)
CD_REAL_420=c(CD_REAL_420,CD_real)
SIMU=SIMU+1;if(SIMU%%1000==0)print(SIMU)
}
write.table(CD_REAL_420,'CD_REAL_420.txt',row.names=FALSE,col.names=FALSE)

SIMU=0;CD_REAL_430=numeric()
for(w in 1:10000)
{
y1=PROFILES_G1_430[,w];y2=PROFILES_G2_430[,w]
DIF=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
dif=(x-1)*stap-(y-1)*stap
DIF[x,y]=dif
}}
DIF=sign(DIF)
DEN=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
den=y1[x]*y2[y]*stap^2
DEN[x,y]=den
}}
CD_real=sum(DIF*DEN)/sum(DEN)
CD_REAL_430=c(CD_REAL_430,CD_real)
SIMU=SIMU+1;if(SIMU%%1000==0)print(SIMU)
}
write.table(CD_REAL_430,'CD_REAL_430.txt',row.names=FALSE,col.names=FALSE)

SIMU=0;CD_REAL_440=numeric()
for(w in 1:10000)
{
y1=PROFILES_G1_440[,w];y2=PROFILES_G2_440[,w]
DIF=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
dif=(x-1)*stap-(y-1)*stap
DIF[x,y]=dif
}}
DIF=sign(DIF)
DEN=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
den=y1[x]*y2[y]*stap^2
DEN[x,y]=den
}}
CD_real=sum(DIF*DEN)/sum(DEN)
CD_REAL_440=c(CD_REAL_440,CD_real)
SIMU=SIMU+1;if(SIMU%%1000==0)print(SIMU)
}
write.table(CD_REAL_440,'CD_REAL_440.txt',row.names=FALSE,col.names=FALSE)

SIMU=0;CD_REAL_450=numeric()
for(w in 1:10000)
{
y1=PROFILES_G1_450[,w];y2=PROFILES_G2_450[,w]
DIF=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
dif=(x-1)*stap-(y-1)*stap
DIF[x,y]=dif
}}
DIF=sign(DIF)
DEN=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
den=y1[x]*y2[y]*stap^2
DEN[x,y]=den
}}
CD_real=sum(DIF*DEN)/sum(DEN)
CD_REAL_450=c(CD_REAL_450,CD_real)
SIMU=SIMU+1;if(SIMU%%1000==0)print(SIMU)
}
write.table(CD_REAL_450,'CD_REAL_450.txt',row.names=FALSE,col.names=FALSE)

SIMU=0;CD_REAL_460=numeric()
for(w in 1:10000)
{
y1=PROFILES_G1_460[,w];y2=PROFILES_G2_460[,w]
DIF=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
dif=(x-1)*stap-(y-1)*stap
DIF[x,y]=dif
}}
DIF=sign(DIF)
DEN=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
den=y1[x]*y2[y]*stap^2
DEN[x,y]=den
}}
CD_real=sum(DIF*DEN)/sum(DEN)
CD_REAL_460=c(CD_REAL_460,CD_real)
SIMU=SIMU+1;if(SIMU%%1000==0)print(SIMU)
}
write.table(CD_REAL_460,'CD_REAL_460.txt',row.names=FALSE,col.names=FALSE)

SIMU=0;CD_REAL_470=numeric()
for(w in 1:10000)
{
y1=PROFILES_G1_470[,w];y2=PROFILES_G2_470[,w]
DIF=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
dif=(x-1)*stap-(y-1)*stap
DIF[x,y]=dif
}}
DIF=sign(DIF)
DEN=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
den=y1[x]*y2[y]*stap^2
DEN[x,y]=den
}}
CD_real=sum(DIF*DEN)/sum(DEN)
CD_REAL_470=c(CD_REAL_470,CD_real)
SIMU=SIMU+1;if(SIMU%%1000==0)print(SIMU)
}
write.table(CD_REAL_470,'CD_REAL_470.txt',row.names=FALSE,col.names=FALSE)

SIMU=0;CD_REAL_480=numeric()
for(w in 1:10000)
{
y1=PROFILES_G1_480[,w];y2=PROFILES_G2_480[,w]
DIF=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
dif=(x-1)*stap-(y-1)*stap
DIF[x,y]=dif
}}
DIF=sign(DIF)
DEN=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
den=y1[x]*y2[y]*stap^2
DEN[x,y]=den
}}
CD_real=sum(DIF*DEN)/sum(DEN)
CD_REAL_480=c(CD_REAL_480,CD_real)
SIMU=SIMU+1;if(SIMU%%1000==0)print(SIMU)
}
write.table(CD_REAL_480,'CD_REAL_480.txt',row.names=FALSE,col.names=FALSE)

SIMU=0;CD_REAL_490=numeric()
for(w in 1:10000)
{
y1=PROFILES_G1_490[,w];y2=PROFILES_G2_490[,w]
DIF=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
dif=(x-1)*stap-(y-1)*stap
DIF[x,y]=dif
}}
DIF=sign(DIF)
DEN=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
den=y1[x]*y2[y]*stap^2
DEN[x,y]=den
}}
CD_real=sum(DIF*DEN)/sum(DEN)
CD_REAL_490=c(CD_REAL_490,CD_real)
SIMU=SIMU+1;if(SIMU%%1000==0)print(SIMU)
}
write.table(CD_REAL_490,'CD_REAL_490.txt',row.names=FALSE,col.names=FALSE)

SIMU=0;CD_REAL_500=numeric()
for(w in 1:10000)
{
y1=PROFILES_G1_500[,w];y2=PROFILES_G2_500[,w]
DIF=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
dif=(x-1)*stap-(y-1)*stap
DIF[x,y]=dif
}}
DIF=sign(DIF)
DEN=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
den=y1[x]*y2[y]*stap^2
DEN[x,y]=den
}}
CD_real=sum(DIF*DEN)/sum(DEN)
CD_REAL_500=c(CD_REAL_500,CD_real)
SIMU=SIMU+1;if(SIMU%%1000==0)print(SIMU)
}
write.table(CD_REAL_500,'CD_REAL_500.txt',row.names=FALSE,col.names=FALSE)

